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Abstract 

We extend next-to-leading logarithmic threshold and joint resummation for prompt photon 
production to include leading collinear effects. The impact of these effects is assessed for both 
fixed-target and collider kinematics. We find them in general to be small, but noticeable. 



1 Introduction 



(1) 



The perturbative QCD description of many observables measured at colliders is plagued by large 
corrections arising from soft and collinear parton emission, even for fairly generic kinematical 
conditions. For example, near threshold, large logarithmic corrections remain [TJ [2] after can- 
cellation of singular virtual and real gluon contributions, their large size being a result of the 
nearby threshold restricting the real gluons to be soft. In terms of a (Mellin) variable N, in 
terms of which threshold is approached by N — > oo, such large threshold corrections take the 
form (L = IniV), 

2i 
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where the depend in general on the process. Another example [3j 01 El [7J [8] is when an 
identified part Fp of a final state has acquired small transverse momentum by soft recoil (Qt) 
against the remaining, unmeasured part of the final state. Then the perturbative expression for 
the differential cross section with respect to px of Fp takes again the form of Eq. ([T]), but with 
different coefficients au and with L = In 6, b being the impact parameter Fourier conjugate to 
Qt- 

Such large logarithmic corrections can be brought under control by all-order resummation, 
and there exists a large literature demonstrating the viability, where applicable, of threshold, 
recoil as well as their joint resummation, for a wide variety of observables. It is interesting 
to try to extend all-order control to classes of large terms beyond the logarithmic corrections. 
One such new set consists of numerically large constants ( 'V 2 terms" ) originating from the same 
infrared-sensitive regions of those Feynman diagrams that also produce the logarithms [91 [TOl [TT] . 

Another important class of potentially large terms, of soft-collinear origin, can be represented 

as 

"■E^— ■ ( 2 ) 
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Their phenomenological importance was first demonstrated in Ref. |12] in which the leading 
terms j = 2i — 1 were also summed to all orders for Higgs production and Drell-Yan. The 
assessment of these terms was made more meaningful in the context of a complete next-to-next- 
leading order (NNLO) [HI El HU [HI CE7J [THl [19] calculation, and a consistent next-to-next- 
leading logarithmic (NNLL) threshold-resummed result [20J. It is not yet clear how to sum 
next-to- leading terms in ([2]). 

In this paper we examine the impact of the leading terms in ([2]) for a single particle inclusive 
observable, the px spectrum of prompt photons produced in hadronic collisions. We do this in 
the context of both a threshold [HI I22J [231 121 [25] and joint [261 123 EH [29] resummed calculation 
for this spectrum. 

The paper is organized as follows. In section 2 we review briefly the threshold and joint 
resummed prompt photon px distribution. In section 3 we describe and motivate our extension 
to include the leading a^ ln ' N — terms. In section 4 we assess the numerical impact of these 
corrections, and we conclude in section 5. 



2 Threshold and joint resummation for prompt photon produc- 
tion 

We consider the inclusive px distribution of prompt photons produced in hadron-hadron colli- 
sions at center of mass (cm) energy y/S 

h A (PA) + h B (p B ) -> 7(Pc) +X, (3) 
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where Ha. b refer to the two incoming hadrons and X to the unobserved part of the final state. 
The lowest order QCD processes producing the prompt photon at partonic cm energy y/s are 

qijpa) + q(Pb) -> 7(p c ) + g(Pd) , ^ 

g(Pa) + q(Pb) -> tOc) + q(Pd) ■ 

The distance to threshold is customarily measured by the variable 1 — x T , where x T = 4p T /S. 
At the parton level this distance is given by 1 — x\ = 1 — 4p T /s. Below we review the result 
for the joint resummed prompt photon pt distribution. At the end of this section we recall how 
the threshold resummed result may be derived from it. 

The joint resummation formalism for prompt photon production \26\ [27] implements the 
notion that, in the presence of soft QCD radiation with summed transverse momentum Qj. of 
soft recoiling partons, the actual transverse momentum produced by the hard collision is not pt 
but rather p^ = Pt — Qt/2. Stated more precisely: in the context of a refactorization analysis 
|27j one can identify a short-distance process at cm energy Q that produces a prompt photon 
with momentum p ' T in a recoiling frame. One defines accordingly x\ = Ap'\ / Q 2 . The extreme 
situation Qt = %Pt in which all transverse momentum is produced through soft recoil leads to 
a singularity in the short-distance process, which we avoid by imposing an upper limit p, on Qt 
|26| . A recently proposed extension [28] of joint resummation avoids this singularity. 

The joint resummed px distribution of prompt photons in hadronic collisions is written as 



„3 j_.(resum) 4 

PT d(7 AB->*i+x Pt dN 



dpT ~ J 8nS 2 Jc 2m 



fi/A(N, fl F )fj/B(.N, fi F ) 



dx 2 T (x 2 T ) N |M * j( ^ )|2 C^ k \a s (fi),x 
. h _ ™2 



x J fi 2 be ib ' QT exp 



N+l 



4p; 



''(■'I 



(5) 



Let us explain each of the terms on the right hand side of Eq. (JS}. The top line displays the 
moments of standard parton distribution functions, as well as the sum over initial state parton 
flavors. The next line contains the Mellin transform over the partonic scaling variable x\ in the 
recoiling frame, the Born amplitudes, and the N- and ^-independent hard virtual corrections 
summarized in C^~^ k K The second to last line contains the integral over the recoil momentum 
of the soft partons, as well as a kinematic factor linking recoil and threshold effects. The last line 
contains the Sudakov exponentials from initial and final state partons, as well as soft wide-angle 
radiation in combined Mellin- impact parameter space. 

As indicated in the last line of Eq. (j5|). large threshold and recoil logarithms, expressed 
through In N and In b, can be resummed into an exponential form. The perturbative exponential 
moment dependence at next-to-leading logarithmic (NLL) accuracy is given by 

Ef T (N, b, Q, M , fi F ) + Ej T (N, b, Q, fx F ) + F k (N, Q, fx) + G ijk (N, M ) . (6) 

Let us discuss each of these terms in turn. The initial state perturbative exponent reads, in 
integral form 

Ef T (N, 6, Q, Hf) = - 2 -jjx- { M".(kr)) In } - 21niV |f ^ MM) 
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where h,hf are the renormalization and factorization scale, respectively. The function x(iV, b) 
defines the N- and 6-dependent scale of soft gluons to be included in the resummation, and is 
chosen as [30] 

x (N,b) = b+^U, (8) 

where rj is a suitably chosen constant and 

N = Ne rB , b = bQe IE /2, (9) 
with "fE the Euler constant. An older form used in [26] 

x (N,b) = b + N (10) 



generates spurious subleading logarithms in Qt [30]. We postpone elaborating on the integral 
in Eq. (|T|) to the next section. The final state jet exponent reads to NLL accuracy 

F k {N,Q,n) = -^—ff\\) + fl 1 \x,Q, f i), (11) 

Qt s \H) 

where 

A = b a s (fi 2 ) \nN . (12) 
The exponent associated with wide angle soft radiation is 

Gabc(N) ^ g^(\) . (13) 

The functions fj^ and g^l(X) as well as the functions C^-? - * 7 ^ [22l [23] are listed in the 
Appendix. 

A nonperturbative term must be added to the perturbative exponent in Eq. ([6]) in order to 
regularize the limit in which Qt is very small. As in Refs. [261 [27] we take 



E^ P = -^ P b 2 ij = qq,qg. (14) 



The threshold-resummed result can now be derived by simply neglecting Qt in the factor 
(5/[4|p T - Q T /2\ 2 ]) N+1 in Eq. ©. Then the Q T integral sets b to zero everywhere, yielding 
the threshold-resummed result. 



3 Including leading In N/N terms 

The leading terms in Eq. ([2]) originate from both initial and final state radiation, and to resum 
them we will use two different methods upon which we elaborate in this section. There are 
moreover two classes of functions in momentum space at order oP s that generate the leading 
ln 2 ^ 1 N/N terms upon Mellin transformation. In terms of the variable z, < z < 1 which in 
the present case can represent either x\ or x T , one of the two classes is formed by the singular 
plus distributions [ln 2:,_1 (l — z)/(l — 2)]+, the other by the singular but integrable ln 2 - ?_1 (l — z). 
The In N/N contributions from the former can be computed using the methods of [31] and can 
be found e.g. in Ref. [12j . The In N/N contributions from the latter can be generated at any 
order in perturbation theory by a simple replacement in the resummed expression (see below in 
Eqs. (j25l) ). expanding the resulting expression to the desired order, and keeping the leading term 
in Eq. ([2]). Roughly speaking, the replacement is equivalent to exchanging at order j one soft- 
collinear gluon (corresponding to one factor a s In 2 N) for a hard-collinear one (corresponding to 
a factor a s In N/N) 

1 2k-l AT 

a k s ln 2k N ^ a k s 1 ^—^ . (15) 
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This replacement is in fact easily included in the existing threshold resummation formulae. A 
preliminary study for prompt photon production was carried out in Ref. [32j. We employ this 
replacement method in fact for the final state related a(?ln 2fc_1 N/N terms. It was pointed out 
in Refs. [30], [33] that the initial state related a k s ln 2fc_1 N/N terms could be generated in the 
context of joint resummation by extending evolution of parton densities to a soft scale. We will 
use this method as well, for the first time for a one-particle inclusive observable. We now discuss 
the initial and final state In N/N contributions in turn. 



3.1 Initial state 

Our procedure for the initial state follows Refs. [30\ 133] . where the joint resummation for elec- 
troweak or Higgs boson production at mass Q and transverse momentum Qp was given. We 
recall the key points here. The integral form of the initial state NLL exponent ([7|) can be written 

as 



Ef T (N,b,Q,H,HF) = - £ | Ai {a s (k T )) In (j^J + B { (a a (fcr))| 

f Q2 /x 2 dkl 

+ -f { -hxNAiMkr^-Biia^hr))} ■ (16) 



The first term in this expression leads to 

Ef T (N,b,Q,v) = -±—h^(f3) + ^\(3,Q,v) , (17) 
where 

P = b a B (ji)hx( X ) ■ (18) 

We recall that the x depends on N and b through Eq. (jSJ). The functions hf 1 ' 1 ^ are listed in the 
Appendix. 

The second term represents flavor-conserving evolution to NLL accuracy (the integrand con- 
sists of the In iV and constant terms for the anomalous dimension matrix 7^ (A) for j = i) from 
the hard scale fip to the soft scale Q/x- One now performs the replacement [301 GS] 

- Ma s ) In (N) - Bi(a t ) — ► 7i/l (N)(a s ) , (19) 

that includes the leading, flavor-diagonal InN/N effects generated by the &t integral (the 1/N 
part of 7j/j combines with the In A" terms). In fact, one may go further and include the off- 
diagonal contributions via the replacement 



5 ig exp 



-A ( g ] InN-B^ 
2^b~ 



f g / H {N, fi F ) — > £ ik (N, Q/x, Mf) fk/ H (N, fi F ) . (20) 



where s{(3) = ln(l — 2/3) plus NLL corrections. As a result, we can replace in Eq. ^ the 
combination 

fi/A^F,N)f j/B ^ F ,N) exp [Ef T (N, b, Q, fi, n F ) + Ej T (N, b, Q, fi, fji F )] (21) 

by 

C i/A (Q, b, N) C j/B (Q, b, N) exp [Ef T (N, b, Q) + Ej T (N, b, M , Q)] (22) 

where 

C i/H {Q, b,N) = J2 £ik (A, Q/x, »f) f k / H (N, hf) ■ (23) 

k 

The matrix £ implements evolution from scale fj, F to scale Q/x, an d is normalized to be the 
unit matrix if these two scales are equal. Note that the dependence on \ip cancels among the 
factors in Eq. (|23[) . 



5 



3.2 Final state 

Leading In N/N effects arising from final state radiation can be derived from the jet functions 
[U H2] that enter threshold or joint resummed expressions for observables having final state 
partons at lowest order. The integral form of the final state exponent F k in Eq. © reads 



1 , z N - 1 - 1 

dz 

1 - z 



^ T A k (a s (q 2 ))+B k (a s ((l-z)Q 2 )) 

(l-z) 2 Q 2 Q 



(24) 



To include leading In N/N dependence in F k (N,Q,fi) we make the replacement [1~2 | 120 1 [52] 



1 - z 



V N-1 



1 



PiZ 



JV-1 



1 



4 1] + o( W2 ) , 



1 - z 

where p q = l,p g = 2. The extra terms can be cast in a more convenient form. Using 

t N ~ x - 1 - (z N - 1) 



JV-l 



1 - z 



and the replacement (accurate to NLL) 



(25) 



(26) 



/- 1 - 1 



1 — jer =■ 

N 



one finds 



F k (N,Q,») 



a s (ix) 



fk\ X ) + fk\ X > Q> + /a(A, « s ) + 0(a,(a s In -/V) n ) , 



(27) 



(28) 



where the extra terms ff k that include the leading In N/N terms due to final state radiation 
read 



A {1) 

J a 



2vr6 



cxp 



A 



34 (1) 
Jg ~ 2irb 



cxp 



a s b 

A 

a s b 



[ln(l - 2A) - ln(l - A)] , 
[ln(l - 2A) - ln(l - A)] . 



(29) 
(30) 



There is no leading In N/N contribution arising from wide angle soft radiation. 

As a result, we finally arrive at the following equation for the joint resummed prompt photon 
hadroproduction px spectrum in which leading soft-collinear effects are included: 



„3 jJresum) 

P T acr AB ^ 



dp T 



i 

Pt 
8ttS 2 



1 ^2 (4)^ \^ 0L>1 C^ k \a s ^),x 2 T ) 



S 



4|pt - Qt/2| 2 



7V+1 



>X i/A (Q, b, N) C j/B (Q, b, N) exp [Ef T (N, b, fi, Q) + Ej T (N, b, fj,, Q)] 



x cxp 



a 8 (p) 



fi° ] W + Q, fA) + f' k (X, a.) + ff gi(A) 



(31) 



As before, the corresponding threshold result may be obtained by neglecting — Qr/2 in the last 
factor on the second line. 
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4 Results 



Here we study numerically the inclusion of the In N/N terms for the case of prompt photon 
production for two kinematic conditions: those of pp collisions at the Tevatron at y/S = 1.96 
TeV (HU [36] , and those of the pN collisions in the E706 [37] fixed target experiment with -Ebeam = 
530 GeV. Our main aim is to assess the effect of such terms in relevant kinematic conditions, 
rather than provide optimized and realistic theoretical calculations for comparison with data (see 
Ref. [38] for recent study). For instance, we do not include contributions from fragmentation 
processes, which have recently been addressed in Ref. |39j and shown to be significant. Our 
assessments mainly consist of comparing the same calculation with and without In N/N terms. 

Our default choices for various input parameters are as follows. We use the GRV parton 
density set [3D] , corresponding to a s (Mz) = 0.114, with the evolution code of Ref. [H], chang- 
ing flavor number at fj, = m c (1.4 GeV) and mf, (4.5 GeV). We choose the factorization and 
renormalization scale equal to pt, and the non-perturbative parameter g-^p in Eq. (jl4H equal to 
1 GeV 2 . For the parameter x we use the expression in Eq. ([8]) , following [30] , with r\ = 1/4 0. 
For our joint-resummed results, we chose for Tevatron (E706) kinematics the cut-off p, in Eq. ([5]) 
equal to 15 (5) GeV. Regarding logarithmic accuracy, and unless otherwise stated we refer to 
LL when using only in Eq. (fT7|) . in Eq. (|28l) . and CK^"^) = 1 for the processes in (j4]); 
we refer to NLL when also including and fjp and the virtual corrections in (|44p . For the 
evolution from scale \ip to Q/x i n Eq. (|20j) we use the full NLO anomalous dimension in all 
cases. 

Starting with Tevatron kinematics we compare in Figs. [T][3] results at LL and NLL accuracy, 
with and without the leading In N/N contribution for joint resummation. For clarity we have 



1.2 - 

1 ^r^—^^ = 

0.8 - """"" 

0.6 - 

4 I — I 1 1 1 1 1 1 1 

20 40 60 80 100 120 140 160 180 20 40 60 80 100 120 140 160 180 

p T [GeV] p T [GeV] 

Figure 1: In N/N contributions for Tevatron kinematics. Left pane: LL without In N/N (a, solid), NLL 
without InN/N (b, dashed), NLL with InN/N (c, short-dashed). Right pane: ratio of NLL to LL (solid), 
ratio of NLL with In N/N to NLL without (dashed). 

here included the constant corrections in (|44|) also for the LL case. Fig. [T] shows that the effect 
of the leading In N/N is appreciable when compared to the effect of passing from LL to NLL, 
the latter difference being almost negligible. Inclusion of In N/N effects leads to noticeable 
suppression for most of the pt range, and to enhancement at very small and very large pr- To 
better understand the origin of these In N/N suppressed contributions, we examine in Figs. [2] 
and [3] for each channel in Q the contributions from the initial and final state. We plot these 
contributions for the LL cross sections only to facilitate interpretation. To help understand the 
results, we can expand the perturbative exponent in Eq. ([6]) to lowest order in a s , keeping only 

1 Choosing rj = 1 does not substantially modify results, but choosing the form in Eq. (|10[) . which generates 
spurious subleading recoil logs [30] , does lead to significant changes at larger pr ■ 
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Figure 2: hi N/N effects for qq channel at LL, Tevatron kinematics. Ratio to LL without \nN/N of 
initial state (solid) and final state (dashed) effects, and both (short-dashed). 
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Pt [GcV] 

Figure 3: hi N/N effects for qg channel at LL, Tevatron kinematics. Labels as in Fig. [2] 

the In 2 N and In N/N terms 

q q : : Si ^ N (2^ - |A<") + S ™ (sMf) - |A<'>) (32) 
qg : : S (4.) + AW _ .4)) + S (4) + 34') - 14)) (33) 

The expressions suggest that the initial state In N/N terms enhance the cross section for the 
qq and in particular the qg channels, while the final state In N/N terms suppress it, again by 
an amount that depends on the channel. The net result turns out to be suppression in the 
former channel and enhancement in the latter. These qualitative aspects are indeed borne 
out if we use the same method to compute initial state In N/N effects as we did for the final 
state in section 13.21 q In the present case however, the net In N/N effect in both channels 
is suppression, indicating that the non-diagonal terms in the evolution matrix give a sizeable 
negative contribution. Note that for Tevatron kinematics, when combining channels, the qg 
channel dominates at low pr, because the required momentum fractions are not too large. At 
large pr, where parton momentum fractions are larger, the valence-quark dominated qq channel 
takes over. 

Turning to E706 kinematics we perform the same studies as we did for the Tevatron. The 
results are shown in Figs. 0H61 We observe an overall enhancement due to the In N/N effects, 
somewhat smaller than the change from LL to NLL. Both effects are more pronounced than for 
the Tevatron. This is due both to a larger value of a s as well as being closer to threshold in 

2 The net result in the qq is actually still enhancement, because the contribution of the f' q ^ g functions in (|29|l . 
()30|) is very small. 
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this fixed target kinematical regime. Examining the effects per channel in Figs. [5] and EJ we now 
see a noticeable enhancement from the initial state In N/N effects in the qq channel, but still 
suppression in the qg channel. Clearly the non-diagonal terms in the evolution matrix play a 
significant role for the E706 case as well. 




1 r 



4 5 6 7 8 9 10 11 

Pt [GeV] 



Figure 5: In N/N effects for qq channel at LL, E706 kinematics. Labels as in Fig. O 
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Figure 6: In N/N effects for qg channel at LL, E706 kinematics. Labels as in Fig. [2 

Next, we examine the differences between threshold and joint resummation. In Fig. [7] we 
compare resummed results directly by showing the ratios with respect to the joint-resummed 
Pt distribution without \nN/N terms. We see for Tevatron kinematics that the threshold 
resummed dominates the joint resummed at large pt, while at low px the converse is true. For 
the E706 case the threshold resummed results are entirely below the joint-resummed ones. The 
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Figure 7: Comparison of joint resummation and threshold resummation effects, ratios to NLL without 
In N/N for Tevatron (left pane) and E706 (right pane). 



threshold resummed curves are shown separately in Fig. which is analogous to the rightmost 
panels in Figs. [T] and [H For Tevatron kinematics the inclusion of In N/N terms in threshold 
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p T [GeVj p T [GeV] 

Figure 8: \nN/N effects in threshold resummation, for Tevatron (left pane) and E706 (right pane). 
Labels as in Fig. [1] right pane. 

resummation leads, as for joint resummation, from suppression at small px to enhancement at 
larger px, but more noticeably. For E706 kinematics, different from the joint resummation case, 
the enhancement at small pt turns to suppression just below px = 6 GeV. The cross section even 
becomes negative beyond 6.5 GeV, which is due to the fact that the nearness of the threshold 
drives the scale Q/x in Eq. (|2U|) effectively below the starting scale of the PDF evolution. 



5 Conclusions 

We have examined the effects of including terms of the form 

■ ( 34 ) 

3 

in joint-resummed and threshold-resummed prompt photon px distributions at both collider 
and fixed target kinematics, at leading accuracy (j = i). The complete structure of subleading 
terms of the form (1341) is still unknown. Note that we have not considered the fragmentation 
component of the prompt photon production cross section in our analysis H. 

3 To do so would require inclusion of more partonic subprocesses, each containing a sum over color structures for 
the wide-angle soft radiation component, as well as photon fragmentation functions [39]. Presumably, soft-collinear 
effects for the fragmentation component of prompt photon production could be included in a way analogous to 
what we did in the present paper for the initial state: via adjustment of the resummed part, and evolution of the 
fragmentation functions. 
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To the extent that terms of the form (|34p arise from initial state radiation effects, we used 
the method of Refs. [301 l33l| to include them, now in a single-particle inclusive cross section. 
Those arising from final state emission we included by extending the jet function to leading 
In N/N accuracy. Numerically we found the combined In N/N terms to be comparable to NLL 
corrections, and dependent on kinematics either enhancing or suppressing. The final state 
In N/N contributions were particularly small, while in the initial state the effects of non-leading 
1/N effects are appreciable, depending again on channel and kinematics. The flavour non- 
diagonal terms in the evolution matrix were found to be numerically significant, and the main 
source of discrepancy with expectations based on simple approximations. We conclude that, 
because the effects, though small, are non- negligible, understanding the structure of In N/N 
terms better is a worthwhile pursuit. 
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A Appendix 

Here we list the exponents used in section The initial state exponents (I17p involve 



A 



(i) 



[2/3 + ln(l-2/3)] 



(35) 



a2\ 

2vr6g 



I ln 2 (1 _ 2/3) + ^±MlzM 

2 y yj 1-2(3 



+ 



B (i) 
vr6n 



ln(l - 2(3) 



2vr6 n 



A™ In 



Here 



where C q = Cf and C g = Ca- Also 



A 



(2) 



7T&0 



1-2(3 



+ ln(l - 2(3) 



C A \ 



67_^ 
18 6 
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T R N F 



(36) 



(37) 



The final state exponents © involve the functions 



(i) 



2vr6 A 



-irbn 



[(1 - 2A) ln(l - 2A) - 2(1 - A) ln(l - A)] 



(38) 



(39) 



f(2) 
J a 



^Ji[ln(l - 2A) - 21n(l - A) + \ ln 2 (l - 2A) - ln 2 (l - A)] 



+ ^ln(l-A) 

7T0 



A 



(2) 



[ln(l - A) - ln(l - 2A)] - ^ [21n(l - A) - ln(l - 2A)] 



4 (1) 



Q 2 



+ ^-[21n(l-A)-ln(l-2A)]ln^ (40) 
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The wide-angle soft radiation exponents (|13p are 

9 ^ {X) = "^ ln(1 " 2A)ln2 ' = "S ln(1 " 2A)ln2 (41) 



In these equations 



IICa-^T r N f nC 2 A -10C A T R N F -6C F T R N F 

12^ ' h = 24^ • (42) 



where Tr = 1/2. These expressions are obtained by expanding the perturbative functions 
A a (a s ), B d (a s ) and D ab ^ dl in powers of a s , 

A a (a s ) = ^AM + (^) 2 AV+0(al) (43) 



7T V 7T 



and so on. 

Finally, the explicit forms of C^~^ k ) [22l [23] are 



- C\ 

IT 



l -(2C F - C A )\n2 + ^K-K q + 2C(2)(2C F - X - 



5 1 3 2« 2 

+ -(2C F - -C A ) In 2 2 + -CH- In 2) - vr6 In -if 



(45) 



7T 



1 1 < A ('2"l 

— (C F - 2C A ) In 2 - -K q + ^(2C F + 19C A ) (46) 



1 3 2r> 
+ -C F In 2 2 + -(Cf + 7r6o)(- In 2) - 7r6 In 

2 4 /i^ 



(47) 



where 



We note that there is no factorization scale dependence in and the coefficient functions in 
Eq. (jUj) because of complete evolution from scale /zp to Q/x in Eqs. (|2Tj) . (|22|) and ([23]) . 
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